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Abstract 

We point out another important production channel of a light pseudoscalar Higgs boson a\ via 
the decays of neutralinos, including X23 ^ Xi a i> i n the framework of the NMSSM. We scan the 
whole parameter space using the most up-to-date version of NMHDECAY and search for regions 
where B(x2 3 ~^ > 0.5. ^ the gluino and squarks are light enough for copious production of 

SUSY events at the LHC, there would be numerous number of X2 3 111 subsequent decays of gluinos 
and squarks. Therefore, the production rates of a\ via neutralino decays would be more important 
than h — > a\a\ and others. Potentially, the final state is filled with many r leptons, which can 
be reconstructed at the mass of a±. This is a clean, observable, and distinguishable signature for 
NMSSM. 

PACS numbers: 12.60.Fr, 13.85.Rm, 14.80.Cp 
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I. INTRODUCTION 



Supersymmetry is the leading candidate for the physics beyond the standard model (SM). 
It not only solves the gauge hierarchy problem, but also provides a dynamical mechanism for 
electroweak symmetry breaking and a natural candidate for the dark matter. The minimal 
version, the minimal supersymmetric standard model (MSSM), has suffered from what so 
called little hierarchy problem and the \i problem. 

An extension with an extra singlet superfield, known as the next-to-minimal supersym- 
metric standard model (NMSSM) was motivated to provide a natural solution to the /i 
problem. The [i parameter in the term nH u H<i of the superpotential of the MSSM naturally 
has its value at either Mpi anc k or zero (due to a symmetry). However, the radiative elec- 
troweak symmetry breaking conditions require the /i parameter to be of the same order as 
rriz for fine-tuning reasons. Such a conflict was coined as the /i problem In the NMSSM, 
the n term is generated dynamically through the vacuum-expectation- value (VEV), v s , of 
the scalar component of the additional Higgs field S, which is naturally of the order of the 
SUSY breaking scale. Thus, an effective \i parameter of the order of the electroweak scale 
is generated. 

The NMSSM was recently revived because it was shown that it can effectively relieve the 
little hierarchy problem Q. Due to the additional Higgs singlet field and an approximate 
PQ symmetry, the NMSSM naturally has a light pseudoscalar Higgs boson a\. It has been 
shown [2I that, in most parameter space that is natural, the SM-like Higgs boson can decay 
into a pair of light pseudoscalar bosons with a branching ratio larger than 0.7. Thus, the 
branching ratio of the SM-like Higgs boson into bb would be less than 0.3 and so the LEPII 
bound is effectively reduced to around 100 GeV jjj. 

The dominance of h — > a^ai mode for the intermediate Higgs boson has significant impacts 
on the Higgs search strategies 4[]. The most useful channel for intermediate Higgs boson, 
h — > 77 and h — > WW* would be substantially affected because B(h — > 77) lowers by a 
factor of a few. So is the h — > bb in Wh, Zh production. New search modes via h — > a±a\ 
are mandatory. For example, h — > a\a\ — > 4b for m ai > Irn^, via Wh, Zh production with at 
least one charged lepton and 45-tags in the final state has been shown to afford a clean signal 
of high significance and a full Higgs mass reconstruction at the LHC [5L The associated 
Higgs production with gauge bosons or tt pairs was shown to be effective [6| . Similar studies 
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h — > a\a\ — > 4r 



at the Tevatron were also performed [7|. Other possibilities like h — > a±a\ — > 2b2r |8j and 



9] can further enhance the signal, especially when 2m T < m ai < 2mb. In 



the extreme limit of zero mixing with the MSSM pseudoscalar, the singlet-like a% can decay 



into a pair of photons. In this case, h — > c^ai — > 47 



10] . The light pseudoscalar _ai can 
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also be produced in non-Higgs decays. It can be produced in B meson decays 
in T decays fl and other rar e decays Q, and in associated production with chargino pan 
M. In some other contexts, a light pseudoscalar boson can also be frequently produced 
in association with a Higgs boson or a heavy quark pair [nj]. There could also be other 
unconventional decay modes of the SM-like Higgs boson in the NMSSM, e.g., invisible decay 
into neutralinos [3]. A recent summary can be found in Ref. 1q| . 

In this note we point out another important production channel of a light pseudoscalar 
Higgs boson a\ via the decays of neutralinos, including x® 3 — ¥ Xi a i? i n the framework of the 
NMSSM. This is potentially much more important than from the decay of the Higgs boson 
or from the associated production. In particular, if the gluino and squarks are light enough 
for copious production of SUSY events at the LHC, there would be numerous number of 
X2 3 m subsequent decays of gluinos and squarks. Therefore, the production rates of a\ via 
neutralino decays would be much more important than h — > a^i and others. There is also 
the possibility of X23 ~^ Xih followed by h — > aia\. It was argued in Refs. {2] that the 
mass m ai in most favorable parameter space is lighter than 2m&, and thus a x — > t + t~ is the 
most frequent. In this case, SUSY events would be filled with many r leptons plus missing 
energies, with the corresponding r + r~ reconstructed at the a\ mass. 

We scan the whole parameter space using the most up-to-date version of NMHDECAY 
[20I and search for regions where B(x%% ~^ Xi a i) > 0-5- We show characteristics of this 
region of parameter space. 



II. TWO BODY AND THREE BODY DECAYS OF NEUTRALINO 

The superpotential of the NMSSM is given by 

W = h u QH u U c -h d QH d D c -h e LH d E c + XS H u H d + \s\ (1) 

where Q, L, H u , H d , U c , D c , E c , and S are the doublet quark and lepton, doublet up- 
type Higgs and down-type Higgs, singlet up-quark and down-quark, and the singlet scalar 
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superfields, respectively. 

The Higgs sector of the NMSSM consists of the usual two Higgs doublets H u and Hd 
and an extra Higgs singlet S. The extra singlet field is allowed to couple only to the Higgs 
doublets of the model, the supersymmetrization of which is that the singlet field only couples 
to the Higgsino doublets. Consequently, the couplings of the singlet S to gauge bosons and 
fermions will only be manifest via their mixing with the doublet Higgs fields. After the 
Higgs fields take on the VEV's and rotating away the Goldstone modes, we are left with a 
pair of charged Higgs bosons, 3 real scalar fields, and 2 pseudoscalar fields. In particular, 
the mass matrix for the two pseudoscalar Higgs bosons P\ and P2 is 

^pseudo = l (Pi P2)M 2 p Pl I (2) 
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with 



Mlu = M\ , 



M 2 P12 = M 2 P21 = ^cot(3 s (M^ sin 2/3- 3A/w 2 ) , 

M P22 = - sin 2/5 cot 2 f3 s (m\ sin 2/5 + 3A/w 2 ) 
3 



M\ = (V2A X + K0 a ) , (4) 



where 

M 2 . = 

sin 2/5 

and tan/5 = v u /vd and tan/5,, = v s /v and v 2 = f 2 + fj. Here Pi is the pseudoscalar in MSSM 
while P 2 comes from the singlet S and from the effects of rotating away the Goldstone modes. 
The pseudoscalar fields are further rotated to the diagonal basis (Ai, A 2 ) through a mixing 



angle (211 ] : 




(5) 



' A 2 \ / cos 8a sin 9 a 
, Ai J y — sin 9 a cos 9a 

where the masses of Ai are arranged such that < m^ 2 . At tree-level the mixing angle 
is given by 

M 2 P12 1 M 2 sin2/5-3AKw 2 
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We also use a\ to denote the A±. The a± is mainly the singlet when 9 a is small. The 
couplings of a\ to fermions scale with sin#A while a\ can have large couplings to Higgsinos 
and Higgs bosons via the XSH u H d term of the superpotential. 



We use the publicly available code, NMHDECAY [20(] , to generate parameter space points. 
Currently, the code has imposed a number of experimental constraints, including the radia- 
tive b — > S7 decay, the Bd and B s mixing parameters, B s — > fi + fi~ decay, B + — > t + u t 
decay, and the relic density of the lightest neutralino. These experimental constraints can 
be turned on or off. The parameter space points presented in this section satisfy all the 
above constraints. Before we show the decay branching ratios of the second lightest neu- 
tralino, we would like to give the vertex factor that we are considering. The vertex factor 

ai-X°-X? is S iven b y 



£ = ai x°i X°j {§ sin 9 A (N j2 - N n tan0 w ) (N i4 cp - N i3Sf3 ) 



X ■ a (AT AT \ AT J XN j3 N H -KN i5 N j5 \ 
--j= sin 6 A [N j3 cp + N jA sp) N i5 + — J - -= J — cos 6 A + {% <-> j) j . (7) 

The other competing channels for %\ — > Xi + X include X = h\, Z, and 3-body decays via 
off-shell particles. The detailed formulas will be shown in a future publication 22j. The 
3-body decays are suppressed as long as at least one of the 2-body modes are open. We show 
in Fig. CTJ the branching ratios of X2 ~ ¥ Xi + ( a ii h\,Z,S — body) for two sets of parameter 
space points, with varying k and A K respectively. It is clear that B(x2 ~ * Xi + a i) dominates 
in these two sets of points, which satisfy all experimental constraints and relic density of the 
LSP. The parameters k and A K can be kept small by the approximate PQ and R symmetries, 
which guarantee the lightness of a\. 

Experimentally, the decay X2 ~~ > Xi + a i gives very interesting signatures. In typical 
SUSY events, a lot of gluinos or squarks are produced, which subsequently decay into the 
third or second lightest neutralinos, instead of directly decaying into the lightest one. Thus, 
there are numerous second lightest or third lightest neutralinos, which then decays into 
the lightest neutralino and the light pseudoscalar a%. The a% then decays into either bb or 
t + t~ depending on its mass. Therefore, the final state will be filled with many r leptons, 
which can be carefully reconstructed at the mass of a\. This is a clean, observable, and 
distinguishable signature for NMSSM. * 



In gauge-mediated models, excessive production of r leptons is often the signature. However, it can be 
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FIG. 1: Decay branching ratios for X2 ~ * Xi + ( a i> h±,Z,3 — body) with (a) A = 0.15, Ay 



-500 



GeV, A K = -0.1 GeV, tan/3 = 10, ji = -150 GeV, M x = 100 GeV, and M 2 = 200 GeV versus «; 
(b) A = 0.2, A x = -500 GeV, k = 0.03, tan/3 = 5, fj, = -150 GeV, Mi = 100 GeV, and M 2 = 200 
GeV versus A K . 

III. DECAYS OF SQUARK INTO NEUTRALINOS 



Here we show the relative branching ratios of a squark decaying into the lightest and 
second lightest neutralinos 

Ql,r -> Ql,r Xi,2 ( 8 ) 

where we simply assume that the Yukawa coupling of the light quark is negligible compared 
with the gauge couplings. The couplings of QL,R~QL,R~Xi are gi ven by 



9u L u L x° 
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N 



9, 



d L d LX ° 



V2 



C2 



N a 



9u R u R x1 
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d R d RX ° 



V2S 

gtw 2 

V2S 



(9) 



We scan the parameter space points with the requirements defined by NMHDECAY, except 
for the relic density of the lightest neutralino, and B(x2 Xi a i) > 0-5 (details in the next 
section). With these parameter space points we calculate the squares of the relative gauge 
couplings: 

\9q LtR q L ^ R x° 2 l 9q LR c- 



s Ql,rX\ 1 



distinguished between the gauge-mediated models and the NMSSM, because the right t pair in NMSSM 
can be reconstructed at m ai but not for the r leptons in GMSB models. 
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FIG. 2: Relative branching ratios of the left-handed or right-handed squark into a quark and 
the neutralin 
X?ai) > 0.5. 



the neutralinos: |a~ „ ~o/o~ „ ~o| 2 . Under the requirements of NMHDEC AY and B(yo 



which can roughly indicate the relative branching ratios without taking into account the 
masses in the final state. We show in Fig. [2] the squares of relative gauge couplings. It 
is clear that the squarks, whether left-handed or right-handed, can frequently decay into 
the second lightest neutralino, instead of just directly into the lightest one. Therefore, it is 
consistent with what we have pointed out in the Introduction that there are numerous second 
or even the third lightest neutralinos via production of squarks or gluinos in SUSY events. 
They will then decay into the lightest neutralino and the light pseudoscalar Higgs boson a\. 
We will map the regions in parameter space that the branching ratio B(x2 3 — > Xi a i) > 0-5 
in the next section. 
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IV. SCAN OF NMSSM PARAMETER SPACE 



Here we scan for the parameter space using the most up-to-date version of NMHDECAY 
J 

201 ] in the following ranges of parameters 

A: - 0.7, A\ : -1000 - 1000 GeV, 

k: -0.7 - 0.7, A K : -10 - 10 GeV, 

tan (3: 1 - 40, fj, : -500 - 500 GeV, 

Mi: - 1000 GeV, M 2 : - 1000 GeV, (10) 

where Mi and M 2 are the bino and wino mass parameter, respectively. We ran for a total 
of 10 million random points in the parameter space. The successful points have to pass the 
criteria of NMHDECAY, including LEPII bounds, b — > sy, Bd and B s mixing, B s — > yU 4 ' [i~ 
and B + — > t + u t , but not the LSP relic density. From the pool of successful points we then 
calculate the branching ratio of B(x\ ~ * Xi a x) an d pick those with the branching ratio larger 
than 0.1. These points are shown in Fig. [3j It is easy to see those points with the branching 
ratio larger than 0.5 in the figure. We also show in Table [I] the reduction of the number of 
points under the requirements of NMHDECAY, and further under B{x% — ► Xi + a i) > 0-1 
and 0.5, respectively. 



V. CONCLUSIONS 

In this note we have scanned the parameter space of the NMSSM with the help of the 
code NMHDECAY. In a sizable fraction of the parameter space points, shown in Table [H 
the branching ratio of B(x2 ~ ¥ X? + a i) > 0-5- It indicates a potentially more important 

TABLE I: Total number of points used, that after scanned by NMHDECAY, and that after im- 
posing B(x2 ~^ Xi a i) > 0-1 an d 0.5. 

Steps Number of points 

Total used 10,000,000 

Scanned after NMHDECAY 41,318 

B(x°2 - X?ai) > 0.1 3,260 

B(x° 2 ^ X?ai) > 0.5 2,030 
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20 40 60 80 100 120 



m ai (GeV) 

FIG. 3: Scatter plot in the plane of m ai and B(x2 ~~ ^ X? a i)- The points pass the criteria of 
NMHDECAY and the branching ratio larger than 0.1. 

production channel of the light pseudoscalar Higgs boson a\ via the decays of neutralinos, 
including — > Xi a i- ft could be much more important than from the decay of the Higgs 
boson or from the associated production. In particular, if the gluino and squarks are light 
enough for copious production of SUSY events at the LHC, there would be numerous number 
of X2 3 in subsequent decays of gluinos and squarks. The a>\ then decays into either bb or 
t + t~ depending on its mass. Therefore, the final state will be filled with many r leptons, 
which can be carefully reconstructed at the mass of a\. This is a clean, observable, and 
distinguishable signature for NMSSM. There is also the possibility of X2 3 ~ * Xih followed 
by h — > aicti, which again gives rise to multi-r-lepton final states. 
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